GLUON PROPAGATORS AND CONFINEMENT 



A.NAKAMURA 
Faculty of Education, Yamagata University 
Koshirakawa, Yamagata990, Japan 
I -in nl n il. i l.< ni .1 1. 1 ; .mi. rj.ii .1-11 .i' || . 

in collaboration with 

H.Aiso, M.Fukuda, T.Iwamiya, T.Nakamura and M.Yoshida 

Computational Sciences Division, 
National Aerospace Laboratory, chofu-city, Tokyo 182, Japan 

ABSTRACT 

We present SU(3) gluon propagators calculated on 48 x 48 x 48 x N t lattices 
at [3 = 6.8 where Nt = 64 (corresponding the confinement phase) and Nt = 16 
(deconfinement) with the bare gauge parameter, a, set to be 0.1. In order to 
avoid Gribov copies, we employ the stochastic gauge fixing algorithm. Gluon 
propagators show quite different behavior from those of massless gauge fields: 
(1) In the confinement phase, G(t) shows massless behavior at small and large 
t, while around 5 < t < 15 it behaves as massive particle, and (2) effective 
mass observed in G(z) becomes larger as z increases. (3) In the deconfinement 
phase, G(z) shows also massive behavior but effective mass is less than in the 
confinement case. In all cases, slope masses are increasing functions of t or z, 
which can not be understood as addtional physical poles. 

1. Introduction 

The title of this conference, "confinement" , is needless to say the most interesting 
and challenging topics in hadron physics. To clarify the phenomena, many approaches 
have been tried but to our knowledge no confirmative study of gluon propagators was 
reported. 

There are two reasons why the non-perturbative study of gluon propagators is im- 
portant to understand the confinement: (f) gluons control the inter-quark dynamics 
and we expect that their propagators show a peculiar behavior at small momentum 
region. (2) the confinement of quarks may be understood as a linear force among 
them but the understanding of the gluon confinement is far from the satisfactory 
stage. Therefore it gives a new insight for the confinement phenomena if we can 
measure directly gluon propagators at the confinement and deconfinement phases. 

Non-perturbative study of gluon propagators is, however, not easy. Analytical 
methods have limited power and need some ansatz. Lattice study is also not simple 
task since (f ) it requires the gauge fixing and Gribov copies may distort propagators, 
and (2) we must measure short and long range regions to study their features and 
therefore work on a large lattice at high f3. 

Here we report our numerical study of gluon propagators in the pure SU(3) gauge 



theory. Propagators of massless gauge fields must have the form, 



GM = < A^p)A v {-p) > 

= 4(v-(i-«)^) a) 

where a is a gauge parameter. We fourier-transform the p sector into t, 

G^( Px ,0,0;t) (2) 
The spatial momentum is chosen as 

p= (^,0,0) = (^,0,0) (3) 

where N x is the lattice size along x direction. 

We call G xx as a longitudinal propagators and G yy and G zz as transverse ones. The 
transverse propagators, Gt, are independent of the gauge parameter a and behave 
as having a simple pole, i.e., exp(— p x t) } if they have the form of Eq. (1). Due to the 
periodic boundary condition, Gj behaves on the lattice, 

G T (t) = Ccosh( Px (t-N t /2)) (4) 

This is our reference form to study propagators, and we shall discuss deviations from 
it. 



2. Gribov ambiguity and its study on compact U(l) 

In 1972, Gribov pointed out that the gauge fixing in non-Abelian gauge theories 
is not unique, and argued that the proper choice of a configuration among copies is 
important for the study of infra-red behavior of the theory. 1 Singer proved rigorously 
the ambiguity problem for SU(N) on S d , and later Killingback found that the same 
mathematical argument holds also for SU(N) and U(l) on T d , where S d and T d are 
d-dimensional sphere and torus. 2 ' 3 

Gribov copies were first observed on a lattice in Ref.(4) for the compact U(l), 
and later on many evidences and analyses for SU(2), SU(3) and U(l) have been 
reported. 5 ' 6 ' 7,8,9 

The compact U(l) model provides us a good place to test algorithms to measure 
gauge propagators: We know that at f3 > 1.0 they must behave as those of photon, 
i.e., Eq. (4). And it was found that the copies distort the propagators, i.e., if one 
simply fixes the gauge with Wilson-Mandula algorithm, 10 ' 11 the propagators do not 
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Fig. 1. Transvers "photon" propagator including effects of Gribov copies. [3 = 1.1 and the lattice 
size is 12 x 8 x 8 x 24. Dotted line is a free massless propagator. Data are taken from Ref.4 



behave in a correct manner. 4 ' 9 See Fig.l. 

3. Stochastic gauge fixing 

Nowaday gauge fixing is a standard tool on the lattice, but it is hard to control 
Gribov copies. One would say that the following definition gives a unique and good 
gauge condition, 12 ' 15 ' 16 



J2 ReTxU^ 



Max. 



(5) 



But numerically it is difficult to ensure the above condition is fulfilled especially on 
large lattices. Several groups have tried to calculate gluon propagators, 12 > 13 > 14 but 
without taking care of the copies. 



We have proposed ' to employ the stochastic gauge fixing by Zwanziger 

dA^ _ 6S 1 
(ir 8A„ a ^ 



19,20 



(6) 



where r stands for the Langevin time, rj is a Gaussian noise term. The second term 
in r.h.s. of Eq. (6) is a force along a gauge trajectory. Usually this force attracts the 
configuration onto a gauge fixed plane where 
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0, 



(7) 



and together with the third noise term, the covariant gauge fixed state with a is 
realized. The algorithm has a remarkable feature, i.e., the second term of Eq.(6) 
brings the configuration such that 
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Fig. 2. Transverse gluon propagator, Gt(^), on 48 x 48 x 48 x 64 at [3 = 6.8 (confinement phase) 



— VA 2 < infi (8) 

dr x 

> out of 0. 

is a region where the Faddeev- Popov operator is positive. Namely the stochastic 
gauge fixing term is attractive (repulsive) inside (outside) of the Gribov region. 

A compact lattice version of Eq. (6) was proposed 17 and correct "photon" prop- 
agators were observed in the Coulomb phase of compact U(f ) model. 



4. Numerical results 



Now we are ready to calculate gluon propagators for SU(3). We use a parallel vec- 
tor computer, the Numerical Wind Tunnel (NWT) at National Aerospace Laboratory, 
Japan, which has the peak performance of 236 GFLOPS and 35 GByte memory. 

The Langevin time step, At, is set to be O.Of. We perform more than 50000 
sweeps for each Monte Carlo run, about f 0% of which is usually discarded. 

In Fig. 2 we show the transverse gluon propagator, Gjit), in the confinement phase, 
together with the free massless one obtained from Eq.(6). Large deviation appears 
as t increases; the gluon seems to become massive at large t. If we would interpret 
the behavior as summing several mass poles, some of the coefficients C in Eq. (4), 
the spectral functions, should be negative. The effective mass seems to vanish again 
at very large t, although there error bars are very large. 

We show in Fig. 3 and Fig. 4 the transverse propagator measured along a spatial 
direction, z, in the confinement and deconhnement phases, respectively. Here effec- 
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Fig. 3. Spatial transverse gluon propagator, Gt{z), on 48 x 48 x 48 x 64 at [3 = 6.8 (confinement 
phase) 
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Fig. 4. Spatial transverse gluon propagator, Gt{z), on 48 x 48 x A. 
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tive masses increase monotonically as z increases. In the deconfinement region, this 
propagator corresponds the screening effect. In the confinement phase, there is no 
such a screening effect, but the behavior of both propagators is similar. The effective 
mass in the confinement is larger than in the deconfinement. 

In conclusion we have observed gluon propagators from short to large range. To 
our knowledge, this is the first measurement of the quantities with taking care of 
Gribov copies. Data are very encouraging, i.e. error bars are under control and the 
results are highly non-trivial; we are now collecting high statistical data. 
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